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ABSTRACT 

We explore possibilities of collapse and star formation in Population III objects ex- 
posed to the external ultraviolet background (UVB) radiation. Assuming spherical 
symmetry, we solve self-consistently radiative transfer of photons, non-equilibrium H2 
chemistry, and gas hydrodynamics. Although the UVB does suppress the formation 
of low mass objects, the negative feedback turns out to be weaker than previously 
suggested. In particular, the cut-off scale of collapse drops significantly below the 
virial temperature Tvir ~ lO^K at weak UV intensities (J21 ^ 10^^), due to both self- 
shielding of the gas and II2 cooling. Clouds above this cut-off tend to contract highly 
dynamically, further promoting self-shielding and II2 formation. For plausible radia- 
tion intensities and spectra, the collapsing gas can cool efficiently to temperatures well 
below lO^K before rotationally supported and the final II2 fraction reaches ~ 10"'^. 

Our results imply that star formation can take place in low mass objects collapsing 
in the UVB. The threshold baryon mass for star formation is ~ 10^ Mq for clouds col- 
lapsing at redshifts z ^ 3, but drops significantly at higher redshifts. In a conventional 
cold dark matter universe, the latter coincides roughly with that of the la density fluc- 
tuations. Objects near and above this threshold can thus constitute 'building blocks' 
of luminous structures, and we discuss their links to dwarf spheroidal/elliptical galax- 
ies and faint blue objects. These results suggest that the UVB can play a key role in 
regulating the star formation history of the Universe. 

Key words: cosmology: theory - diffuse radiation - galaxies: formation - radiative 
transfer - molecular processes 



1 INTRODUCTION 

The existence of an intense ultraviolet background (UVB) 
radiation, inferred from observations of QSO absorption 
spectra (Gunn & Peterson 1965; Bajtlik, Duncan & Ostriker 
1988), is likely to have prominent impacts on galaxy for- 
mation. While its origin could be attributed to QSOs or 
young galaxies appeared at high redshifts (e.g., Couchman 
1985; Miralda-Escude & Ostriker 1990; Sasaki & Takahara 
1994; Fukugita & Kawasaki 1994), the UVB once produced 
photoionizes the intergalactic medium and alters the subse- 
quent growth of cosmic structures. These secondary struc- 
tures emerged in the UVB are further responsible for pro- 
ducing the radiation field, which can in turn affect the devel- 
opment of next generation objects. It is therefore crucial to 
determine quantitatively the consequences of radiative feed- 
back on the formation of early generation (hereafter Popu- 
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lation III) objects, in order to link the build-up of structures 
with the thermal history of the universe. 

It has long been suggested that the formation of sub- 
galactic objects is suppressed via photoionization and heat- 
ing caused by the UVB (Umemura & Ikeuchi 1984; Ikeuchi 
1986; Rees 1986; Bond, Szalay & Silk 1988; Efstathiou 1992; 
Babul & Rees 1992; Murakami & Ikeuchi 1993; Chiba & 
Nath 1994; Zhang, Anninos & Norman 1995; Thoul & Wein- 
berg 1996, among others). Most of previous studies, however, 
assumed for simplicity that the medium is optically thin 
against the external radiation. Several attempts have been 
made recently to take into account radiative transfer of pho- 
tons (e.g., Kepner, Babul & Spergel 1997; Tajiri & Umemura 
1998; Abel, Norman & Madau 1999; Barkana & Loeb 1999; 
Razoumov & Scott 1999; Gnedin 2000b; Susa & Umemura 
2000a; Nakamoto, Umemura & Susa 2001; Ciardi et al. 
2001). In our previous papers (Kitayama & Ikeuchi 2000, 
hereafter Paper I; Kitayama et al. 2000, hereafter Paper II), 
we have studied the evolution of spherical clouds exposed 



© 2001 RAS 



2 T. Kitayama et al. 



to the UVB, solving explicitly radiative transfer of ionizing 
photons and hydrodynamics. We have found that the col- 
lapse of objects with circular velocities below 14 ~ 15 — 40 
km s~^ can be prohibited completely even if self-shielding of 
gas is taken into account. We have also noted that a spheri- 
cal cloud above this threshold undergoes run-away collapse 
in the dark matter potential and does not settle into hydro- 
static equilibrium. For plausible intensities and spectra of 
the UVB, the central region of such a cloud is shown to get 
self-shielded against photoionization before shrinking to the 
rotation barrier. A question still remains as to whether star 
formation can take place efficiently in the cloud collapsing 
under the UVB. 

A key element in this context is molecular hydrogen 
(H2). In order for stars to form, the gas needs to radiate 
energy efficiently and cool down to temperatures well be- 
low T = 10*K. In the gas of primordial composition with- 
out metals, cooling processes at this temperature are almost 
solely dominated by rotational- vibrational line excitation of 
H2 (e.g., Peebles & Dicke 1968; Hirasawa 1969; Matsuda, 
Sato & Takeda 1969). Formation and destruction of H2, how- 
ever, are very sensitive to the presence of a radiation field. 
For instance, H2 molecules are easily dissociated by photons 
in the so-called Lyman-Werner (LW, 11.2-13.6 eV) bands 
(Stecher & Wilhams 1967; Haiman, Rees & Loeb 1997; Cia- 
rdi et al 2000b) or via collisions with ions. On the other 
hand, the gas-phase formation of II2 is promoted by an en- 
hanced ionized fraction in the post-shock gas (Shapiro & 
Kang 1987; Kang & Shapiro 1992) or in weakly photoion- 
ized media (Haiman, Rees & Loeb 1996b). In addition, self- 
shielding taking place during dynamical collapse can also aid 
H2 cooling (Susa & Kitayama 2000). It is thus by no means 
trivial in what circumstances the UVB has positive or neg- 
ative feedback on star formation in Population III objects. 

In this paper, by including H2 chemistry consistently 
with radiative transfer of UV photons, we extend our pre- 
vious analyses (Papers I and II) and study the fate of Pop- 
ulation III objects exposed to the external UVB. Assuming 
spherical symmetry, we solve the radiative transfer equation, 
non-equilibrium chemical reactions including H2 formation 
and destruction, and dynamics of baryonic gas and dark 
matter. Particular attention is paid to the possibilities of 
both collapse and H2 cooling in the presence of UV radi- 
ation fields. We explore the evolution of objects with wide 
ranges of masses, including those with the virial tempera- 
ture Tvir < lO^K (K < 17 km s"^), which have been left 
out in our previous analyses neglecting H2 formation. A sim- 
ilar problem was investigated earlier by Kepner et al. (1997) 
and Haiman, Abel & Rees (2000), assuming that collapsed 
clouds are in approximate hydrostatic equilibrium. As men- 
tioned above and shown explicitly in this paper, the gas 
in collapsed objects is far from static and its physical state 
changes highly dynamically. We also note that our analysis is 
complementary to that of Susa & Umemura (2000a,b), who 
explored the pregalactic pancake collapse in plane-parallel 
geometry and studied the fate of objects far above the Jeans 
scale. We pay particular attention to the development of 
smaller objects near the Jeans scale, which are likely to 
evolve spherically to a first-order approximation. A crucial 
distinction between the earlier works and ours is that the 
gas density rises rapidly along with the spherical collapse. 



in proportion to the cube of scale length, providing greater 
chances of self-shielding and H2 cooling. 

This paper is organized as follows. Section 2 describes 
the numerical model used in the paper. We present the sim- 
ulation results in Section 3 and discuss their implications on 
galaxy formation in Section 4. Finally, Section 5 summarizes 
our conclusions. Throughout the paper, we assume the den- 
sity parameter f2o = 0.3, the cosmological constant Ao = 0.7, 
the Hubble constant /iioo = //o/(100 km s~^ Mpc"^) = 0.7, 
and the baryon density parameter f2b = 0.05. 



2 MODEL 

2.1 Numerical Scheme 

We simulate the evolution of a spherically symmetric cloud 
exposed to the external UVB. The cloud is a mixture of 
baryonic gas and dark matter with the mass ratio of fib : 
f2o — f2b = 1 : 5. The basic equations are summarized in Sec- 
tion 2.1. of Paper I. One modification is that the acceleration 
term due to the cosmological constant +\(,H^r^ is added to 
the right hand side of momentum equations. They are solved 
with the second-order-accurate Lagrangian finite-difference 
scheme described in Bowers & Wilson (1991) and Thoul & 
Weinberg (1995). Shocks are treated with an artificial ten- 
sor viscosity (Tscharnuter & Winkler 1979). The number of 
mass shells is A^b = 300 for baryonic gas and A'^d = 10, 000 
for dark matter. Our code reproduces accurately the simi- 
larity solutions for the adiabatic accretion of collisional gas 
and the pressureless collapse around a point-mass pertur- 
bation (Bertschinger 1985). We have also checked that our 
results are essentially unchanged when the shell numbers are 
doubled or halved. 

At each time-step of the simulation, we solve self- 
consistently the radiative transfer equation, chemical reac- 
tions, and energy equation in the following three steps. Re- 
peating these procedures, the radiation field and the internal 
state of the gas are determined iteratively until the abun- 
dance of each species and the internal energy in each mesh 
converge within an accuracy of 0.1%. 

First, the direction/frequency-dependent radiative 
transfer is worked out for ionizing photons with energies 
hi/ > 13.6 eV, with the method devised by Tajiri & 
Umemura (1998), in which both absorption and emission 
by hydrogen are explicitly taken into account. At each ra- 
dial point, angular integration of the transfer equation is 
done over at least 20 bins in 9 = — n, where 9 is the an- 
gle between the light ray and the radial direction. This is 
achieved by handling 400-1,000 impact parameters for light 
rays. Self-shielding of H2 against the LW photons (11.2-13.6 
eV) is evaluated from the self-shielding function of Draine 
& Beltoldi (1996), considering the direction dependence of 
the H2 column density as in Kepner et al. (1997). The gas 
is assumed to be optically thin against the photons with 
hf < 11.2eV, which only contribute to destructions of H~ 
and HJ. 

Secondly, non-equilibrium chemical reactions are solved 
in each mesh for the species e, H, H^, H~ , H2, HJ. Hydrogen 
molecules form mainly via the creation of the intermediates 
H" and H+ (eqs [1.1]-[1.4] of Shapiro & Kang 1987). We 
neglect helium for the following reasons: (i) it only changes 
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the ionization degree of the gas up to ~ 10% (Osterbrock 
1989; Susa & Umemura 2000a); (ii) heUum Une cooUng is 
unimportant at the typical temperatures (T 5: lO^K) of the 
objects presently studied; and (iii) photoionization heating 
of helium does not alter our main results as discussed in 
Section 3.1. Unless stated explicitly, reaction rates are taken 



from the recent compilation of Galli & Palla (1998). The co- 
efficients for photoionization of H and photodissociation of 
H2 are computed from the results of the radiative transfer 
mentioned above. Photodetachment of H~ and photodisso- 
ciation of are also taken into account using the cross 
sections given in Tegmark et al. (1997) and Stancil (1994) 
for the photons with hv > 0.74 eV and hu > 0.062 eV, re- 
spectively. These two processes, particularly the former, can 
suppress H2 formation by destroying its intermediates. 

Finally, the energy equation is solved including 
UV heating and radiative cooling due to coUisional 
ionization, collisional excitation, recombination, thermal 
bremsstrahlung, Compton scattering with the cosmic mi- 
crowave background (CMB) radiation, and rotational- 
vibrational excitation of H2. The UV heating rate is ob- 
tained directly from the results of the radiative transfer. The 
atomic cooling rates are taken from Fukugita & Kawasaki 
(1994) and molecular cooling rates from Galli & Palla 
(1998). 

For comparison with previous approaches, we also per- 
form simpler calculations in the following two ways. One 
assumes at each time-step an optically thin medium and 
ionization equilibrium among e, H and H+, neglecting H2, 
H~ and Hj. The other solves radiative transfer of ionizing 
photons, but neglects H2, H~ and HJ. The latter corre- 
sponds to the analysis of Paper II. Hereafter we refer to the 
former as the 'optically thin calculation', the latter as the 
'no H2 calculation', and the one described in the previous 
paragraphs as the 'full calculation'. Unless otherwise stated, 
the full calculation is adopted. 



2.2 External UVB 

The external UVB is assumed to be isotropic and to have 
either a power-law or black-body spectrum: 



7° 



(1) 



xp(/ll//feBTaft)-l' 



where fee is the Boltzmann constant. We adopt the spectral 
index a = 1 and the effective temperature T^a = lO^K, 
to mimic the spectra of QSOs and stars, respectively. The 
above spectrum is normalized to give the intensity J21, in 
units of 10~^'^erg s~^ cm~^ str~^ Hz"^, at the Lyman limit 
of hydrogen {hf — 13.6eV). 

Fig. 1 illustrates the model spectra described above with 
J21 = 1. Also shown for reference are those with a — 5 and 
Tefl = 5 X 10*K. In the present analyses, the ratio between 
the fluxes below and above the Lyman limit is a key pa- 
rameter to control the radiative feedback on the gas evolu- 
tion (e.g., Ciardi, Ferrara & Abel 2000a). To help quantita- 
tive comparisons among different spectral shapes, we list in 
Table 1 the photon-number-weighted average intensities of 
these spectra: 

(J) = — X 




Figure 1. Model spectra of the external UVB with a = 1 (solid), 
a = 5 (dotted), T^ff = IQ-'K (short dashed), and T^s = 5 X 10'*K 
(long dashed), all normalized to give J21 = 1. 



Table 1. Average intensities (J) defined by equation (Q) in spe- 
cific energy ranges. All spectra are normalized to give J21 = 1. 



spectrum 



0.74-13.6 oV 11.2-13.6 oV 13.6-10"' eV 



a = 


1 




6.0 


1.1 


0.15 


a = 


5 




1.4 X 10^ 


1.7 


3.0 X 10-2 






10"* K 


1.3 X 10^ 


3.8 


1.1 X 10-2 






5 X 10-* K 


0.39 


1.0 


8.6 X 10-2 



J 1/hv dv 



10 



'^^erg s ^ cm 



str 



Hz" 



(2) 



where the average is taken over the energies relevant for 
destruction of both H2 and H- (0.74-13.6 eV), dissociation 
of only H2 (11.2-13.6 eV), and ionization of H (13.6-10** eV). 
Notice that for a given value of J21 the spectrum with T^s — 
lO^K has considerably less ionizing (> 13.6eV) photons and 
more H2 and H- destroying photons (< 13.6eV) than the 
one with a — 1. 

Observations of the proximity effect in the Lya forest 
suggest J21 = 10='=° '^ at redshift z = 1.7 - 4.1 (Bajtlik et 
al. 1988; Bechtold 1994; Giallongo et al. 1996; Cooke, Es- 
pey & Carswell 1997), but its value is still highly uncertain 
at other redshifts. Recent CMB power spectrum measure- 
ment by MAXIMA has been used to constrain the redshift of 
reionization to be 2 X!, 8 at the 2a confidence level (Schmalz- 
ing, Sommer-Larsen & Gotz 2000). In this paper, we fidu- 
cially fix the onset of the UVB at zjjv — 20 and consider 
at z < zjjv constant J21 of an arbitrary value. To take into 
account uncertainties in the reionization history, we also ex- 
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amine some cases with zuv = 10 and 50, and evolving J21 
of the form: 

{exp[-(z-5)] 5 < z < 20 
1 3 < 2 < 5 (3) 

{^Y 0<z<3, 

where the evolution at z > 5 is roughly consistent with 
recent models of reionization of the universe (Ciardi et al 
2000b; Umemura, Nakamoto & Susa 2000). 



2.3 Initial and boundary conditions 

Initial and boundary conditions are similar to those in Sec- 
tions 2.3. and 2.4. of Paper I. In the following, we give their 
brief summary and state any modifications and additions. 

We start the simulations when the overdensity of a 
cloud is still in the linear regime. The initial overden- 
sity profile for both baryon and dark matter is Si{r) = 
Si{0) sm{kr)/kr, where k is the comoving wave number, and 
the central overdensity 5i(0) is fixed at 0.2. Only the grow- 
ing mode is considered. The outer boundary is taken at the 
first minimum of 5i{r), i.e., kr = 4.4934, within which the 
volume averaged overdensity Si{< r) vanishes. The charac- 
teristic mass of a cloud Mdoud is defined as the baryon mass 
enclosed within the first zero of 5i(r), i.e., kr — tt, which 
is 0.363 of the total baryon mass (Haiman, Thoul & Loeb 
1996a). We define the central collapse redshift Zco and the 
cloud collapse redshift Zc, respectively, as epochs at which 
the inner-most gas shell and the shell enclosing Mdoud would 
collapse to the center in the absence of thermal pressure. 
These two redshifts are related by 1 -I- Zco — 2.7(1 + Zc) 
(strictly speaking, this relation is correct in the Einstein-de 
Sitter universe, but it still gives a reasonable approximation 
in the present case at Zc ^ 1). We will often characterize a 
cloud by its circular velocity K and virial temperature Tvir, 
defined in terms of Zc and Afcioud as 



Vc = 15.9 



V lO^/irooMo 



Ac(gc) 

187r2 
9.09 X 10^ 



100 ' 
1/6 



1/3 



(1 + %) 



1/2 



km s 



(4) 



0.59 

1/3 



A/doudSlo/^^b 



2/3 



187r2 



(l + 2c) K, 



Mr; 



(5) 



where fj, is the mean molecular weight in units of the proton 
mass rup, and Ac(2:c) is the mean overdensity of a virial- 
ized object predicted from the spherical infall model (Pee- 
bles 1980; Kitayama & Suto 1996). Using circular velocity 
or virial temperature to denote the size of an object, we 
can minimize the dependences of our simulation results on 
assumed cosmological parameters. 

The initial abundance of each species is taken to be its 
cosmological post-recombination value (e.g., Galli & Palla 
1998), i.e., Xo = 10"^, X^- = 10-^^ Xm = 10-^ and 
X^2+ ~ 10"^'', where Xi = ni/riH is the fraction of the 
species i with respect to the total number of hydrogen atoms. 
Our results are insensitive to the precise values of initial 



abundances. In the optically thin and no H2 calculations, 
Xe = lO"** is assumed before the gas is exposed to the UVB 
or experiences shocks. 

In this paper, we mainly focus on the evolution of 
the central gas of a cloud until it reaches the rotation ra- 
dius specified by the dimensionless spin parameter at turn- 
around Ata: 



r,ot = 0.03 



0.17 



Ata 

005 



(6) 



where rta is the turn-around radius of the shell, and we adopt 
a median for the spin parameter, Ata = 0.05 (Efstathiou & 
Jones 1979; Barns & Efstathiou 1987; Warren et al. 1992). 
The above radius defines the scale below which the gas is 
supposed to attain rotational support and the approxima- 
tion of spherical symmetry breaks down. The overdensity 
of the gas exceeds ~ 10^ when it reaches this radius. The 
simulation is halted when the innermost gas shell has fallen 
below this radius and is regarded as 'collapsed'. On the other 
hand, a dark matter shell falling to the center is rebounced 
at the radius Trot of the innermost gas shell, to mimic shell 
crossings. 

At the outer boundary, we assume that the gas is in 
pressure equilibrium, i.e., the pressure outside the cloud Pout 
is equal to that of the outermost gas shell. This is different 
from the free boundary condition (Pout = 0) adopted in Pa- 
pers I and II. While the inner structure of a simulated cloud 
is hardly afi'ected by this difference, we adopt the former in 
this paper so that the cloud envelope maintains to follow 
the Hubble expansion throughout its evolution. 



3 RESULTS 

3.1 Significance of radiative transfer and H2 
formation 

The results of the full calculations are first compared with 
those of previous approaches, to illustrate the significance 
of incorporating radiative transfer and H2 formation in our 
analyses. The most noticeable changes are in the evolution 
of a low mass object with Tvir < lO'^K. Fig. ^ shows the time 
evolution of the central part of a cloud with Zco = 10 (zc = 3) 
and Vc = 10 km s~^, corresponding to Tvir = 4 x lO^K. Both 
the dynamics and the thermal state are shown to change 
completely under different treatments of the radiative pro- 
cesses. In the optically thin calculation, the cloud is pho- 
toionized to the center and evaporated promptly by a little 
radiation with J21 — 10~^ and a = 1. In the no H2 calcu- 
lation, the external radiation is attenuated by absorption, 
but contraction of the gas is halted by thermal pressure as 
it lacks coolant at T < lO'^K. In the full calculation, on the 
other hand, the cloud center is kept self-shielded against the 
UVB and H2 molecules in excess of Xu2 ~ 10"'' are pro- 
duced, allowing the gas to collapse and cool efficiently to 
~ lO^K. Details of the H2 formation and destruction pro- 



cesses in the UVB will be discussed in Section i.2. 

For a high mass object well above the scale Tvir = 10*K, 
in contrast, radiative transfer and H2 formation appear to 
have smaller impacts as far as the dynamics is concerned. 
Fig. ^ shows that such a cloud can collapse even in a stronger 
UV field with J21 = 1 and a — 1, regardless of treatments of 
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Figure 2. Time evolution of (a) radius, (b) temperature, (c) electron fraction, and (d) H2 fraction, of the central gas shell of a low mass 
object with Vc = 10 km and z^o = 10 (zc = 3), comparing different treatments of radiative processes. Lines and symbols indicate 
results of the optically thin (dotted line), no H2 (circles), and full (solid line) calculations with J21 = 10~^ and a = 1. 



radiative processes. This is because gravitational potential 
of the cloud is much greater than thermal energy attained by 
photoionization (T ~ lO^K) and atomic cooling is efficient 
enough to enable the collapse. 

Nevertheless, the final thermal state of the gas is still 
modified to a great extent in a high mass collapse, as in- 
dicated by the differences in the temperature and electron 
fraction at the collapse in Fig. |^. To see this point more 
clearly, Fig. ^ compares the internal structure of the same 
objects at the collapse of the central gas. The density profile 
commonly exhibits that the central gas is undergoing run- 
away collapse (panel a), as it cools efficiently and contracts 
nearly at the free-fall rate under the dark matter poten- 
tial. The radiative transfer effects, however, produce a self- 
shielded neutral region at r ^ 2kpc (panel b), whereas Xm 
is simply determined by ionization equilibrium with uniform 
radiation in the optically thin case. The gas well inside this 
neutral region can cool down to ~ lO^K once H2 cooling is 
taken into account (panel c). In the absence of H2 molecules. 



the gas can cool only to ~ IO^'K by atomic cooling. Rapid 
increase of density due to run-away collapse is in fact es- 
sential for self-shielding and H2 cooling, as will be discussed 
throughly in the forthcoming sections. 

To further corroborate the above results. Fig. ^ shows 
more detailed structures in the full calculation. While H2 
molecules are destroyed in the ionized outer envelope, they 
form efficiently in the self-shielded central region (panel a). 
Roughly speaking, the gas becomes optically-thick against 
ionizing (> 13.6eV) photons at A^hi ~ lO^'^ cm~^ and the 
LW (11.2 - 13.6eV) photons at iVH2 ^ lO^"* cm"^, where 
A''hi and A'h2 are respectively the HI and H2 column densi- 
ties measured from the outer boundary. Both of these con- 
ditions are amply satisfied at r ^ 1 kpc (panel b). The UV 
heating rate is reduced significantly in this optically thick 
region (panel c). Photons emitted in the cloud interior con- 
tribute at most ~ 10% to the total heating rate for the a — 1 
spectrum, though their contribution can be much larger for 
a softer spectrum (Tajiri & Umemura 2000) . Cooling to tem- 



© 2001 RAS, MNRAS 000, [|-|l| 



6 T. Kitayama et al. 




Figure 3. Same as Fig.^ except for showing a high mass collapse with Vc = 50 km s ^ in the UVB with J21 = 1 and a = 1. 



peratures below 10'*K is possible once H2 cooling overtakes 
UV heating. At T < 10''K (r < 0.8 kpc), the timescale of 
H2 formation becomes shorter than that of destruction as 
recombination lags behind the temperature decrease (Kang 
& Shapiro 1992; Corbelh, Galh & Palla 1997), and chemi- 
cal reactions are out of equilibrium (panel d). The cooling 
timescale, however, becomes even shorter as temperature 
drops, and Xb_2 is frozen out at ~ 10~^ at the cloud center. 
This abundance is comparable to that achieved commonly in 
a metal-deficient post-shock layer in the absence of radiation 
(Shapiro & Kang 1987; Ferrara 1998; Susa et al. 1998). 

If helium is further incorporated into the present cal- 
culations, the largest impact will be via heating by photons 
with hu > 24.6 eV. For given UV intensity J21, HI number 
density nni, and HI column density A'hi, the heating rate 
will be enhanced by the addition of helium, in the case of 
a = 1, as 



n{U + He) 



1 + 0.19 ( ""g'os ) optically thin limit. 



1.46 



i/'i 



iVHcl/iVHI 



optically thick limit. 



where 7i(H + He) is the total heating rate in units of erg 
s~^ cm~^ of both H and He, and H(H) is that of the pure 
hydrogen gas. These rates are computed as in the Appendix 
of Paper I, which are in good agreement with more rigorous 
treatment. Equation (j^ indicates that the total heating rate 
tends to increase by ~ 20 — 50% once helium with the pri- 
mordial abundance (~ 8% by number) is taken into account. 
This is much smaller than the overall variations in the heat- 
ing and cooling rates in a cloud shown in Fig. |^(c). In fact, 
we have checked that an additional heat input of 50% does 
not considerably alter either the time evolutions in Figs ^ 
and ^ or the internal structures in Figs ^ and ^. Though 
temperature of the diffuse outer envelope can increase by as 
much as a factor of 2 by helium heating (Abel & Haehnelt 
1999), the evolution of the central high-density region, in 
which our main interests reside, is little influenced. Helium 
heating will be of less importance for softer spectra. Hence 
we neglect helium in our analyses. 

In summary, the above results justify our present ap- 
(7) proach and confirm the importance of including radiative 
transfer of UV photons, H2 formation/destruction, and hy- 
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Figure 4. Radial profiles at the collapse of the cloud center, comparing the optically thin (dotted line), no H2 (circles), and full (solid 
line) calculations: (a) total hydrogen density, (b) HI fraction, and (c) temperature. J21 = 1, o = 1, 14 = 50 km and 2cO = 10 (zc = 3) 
are adopted. 



drodynamics, in studying the development of population III 
objects in the UVB. 

3.2 Dynamical and thermal evolution 

The evolution of UV irradiated clouds are studied in more 
detail in the following. Fig. ^ displays how the time evolu- 
tion of a low mass object (Vc = 10 km s~^) depends on the 
radiation intensity. In the absence of the UVB ( J21 — 0), the 
cloud decoupled from the Hubble expansion is first shock- 
heated and then cools radiatively. At J21 = 10~^, in con- 
trast, the gas is photoionized and evaporated even if ra- 
diative transfer is taken into account. Once the gas is pho- 
toionized and heated, H2 molecules are destroyed by a large 
amount. At lower intensity (J21 = 10"*), as was also shown 
in Fig. ^, the cloud can collapse and cool down by H2. At 
even weaker intensity (J21 — lO"'^), the radiation is attenu- 
ated almost completely and the evolution resembles that in 
the no UV case. Note that the H2 fraction increases after the 



onset of the UVB at J21 = lO"'', compared to the no UV 
case. This is indeed the case in which the formation of H2 
via H~ and is promoted by an enhanced ionized fraction 
in a weakly photoionized medium (Haiman et al. 1996b). 

Fig. 1^ illustrates that the shape of the radiation spec- 
trum also infiuences greatly the thermal and dynamical evo- 
lution. Compared to the a — 1 case, the black-body spec- 
trum with Tss = 10'*K is softer and has much less ionizing 
photons with energies above 13.6eV (see Fig. 1 and Table 
1). For the same values of J21 as in Fig. ^ the black-body 
UVB therefore has weaker impacts on the evolution of ra- 
dius, temperature, and the ionization degree. On the other 
hand, H2 molecules are destroyed by a larger amount at the 
onset of the UVB, because a greater number of photons with 
energies below 13.6eV are available to destroy H2 and its in- 
termediates, H~ and H^. Nevertheless, as the gas contracts, 
it becomes optically thick against the LW photons and the 
final H2 fraction of Xh2 ~ 10~^ is still achieved. 

For a high mass object well above the scale Tvir ~ 10*K, 
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timescales of free-fall (dashed), cooling (dot-dashed), H2 destruction (dotted), and H2 formation (solid). 



the collapse is hardly suppressed even under the strong 
UVB, but its thermal evolution remains sensitive to the 
UVB. As shown in Fig. ]A for J21 ^ 1 and a = 1, the gas is 
promptly heated to ~ 10 K and H2 molecules are destroyed 
mainly via collisions with H+. At J21 = 10~^, on the other 
hand, the H2 abundance is enhanced compared to the no UV 
case by an increased degree of ionization, as already men- 
tioned. In either case, the H2 fraction at the collapse again 
exceeds Xh2 ~ 10~^. 

To illustrate more clearly how self-shielding and radia- 
tive cooling operate in a collapsing cloud, we plot in Figs ^ 
and |lo| the temperature and fractions of HI and II2 at the 
cloud center as a function of its density. As the gas contracts 
and density rises, its state changes rather suddenly from the 
hot (T > lO^K) ionized 'HII phase', to the warm (T ~ lO^K) 
neutral 'HI phase', and the cold (T < 10*K) neutral 'H2 
phase'. Such transitions are analogous to those pointed out 
by Kepner et al. (1997), except that they are induced here by 



the dynamical contraction of the gas itself rather than the 
decline of the external UV intensity. The transition from 
the HII phase to the HI one takes place when the cloud 
becomes optically thick against ionizing photons. The gas 
finally enters the H2 phase once the cooling timescale be- 
comes shorter than the other relevant timescales, such as 
those of free-fall and H2 formation/destruction. These fig- 
ures confirm that Xu2 tends to converge to ~ 10~^ after the 
H2 phase is achieved. 

Figs ^ and ^ further indicate that quantitative details 
of the transitions among the above phases depend sensitively 
on the intensity and spectrum of the UVB. To examine this 
more directly, we plot in Fig. |ll| the critical densities for HI 
shielding (HI) and H2 cooling rig (H2), defined as those at 
which the central HI fraction exceeds 0.5 and temperature 
falls below 5000K, respectively. In what follows, we regard 
the gas as in the HII, HI, and H2 phases, if the density is 
in the range nn < nH(HI), rin (HI) < nn < nH(H2), and 
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Figure 6. Same as Fig.^ except for comparing different values of J21 in the case of o = 1. Lines indicate the results with J21 = 10 ^ 
(long dashed), 10~^ (solid), 10~® (short dashed), and (dotted). The results of J21 = 10~^ and closely overlap each other. 



riH > 71h(H2), respectively. Though Fig. |l^ is plotted for 
particular choices of T4 and Zc, basic features are similar in 
other cases. A possible change is that for smaller masses and 
large J21 the collapse is completely prohibited (see Fig. [l2| ) 
and the critical densities simply disappear. 

Fig. [n] shows that tih (HI) increases with J21 for a given 
spectral shape, i.e., the higher density is required to shield 
the stronger ionizing flux. At small J21, 72^(111) falls below 
the turn-around density, indicating that photoionization is 
kept attenuated throughout the cloud evolution. This fea- 
ture is also apparent in Figs I and 0, in which Xui is kept 
high at weak UV intensities. Compared to the a — 1 case, 
the Teff = lO^K spectrum yields lower values of (HI), be- 
cause there are less ionizing (> 13.6 eV) photons for a given 
J21. 

On the other hand, (H2) behaves in a more complex 
manner, depending on the spectral shape. For q = 1, it is not 
a monotonic function of J21, but has a minimum at moder- 
ate intensities J21 ~ 10"'^ — 10~^. Formation of H2 molecules 
is enhanced maximally at these intensities. At weaker inten- 



sities J21 10~* the result almost coincides with that of the 
no UV case, whereas at stronger intensities J21 10~^ H2 
formation is suppressed to higher densities. For Teff = 10*K, 
in contrast, a greater fraction of photons below 13.6 eV is 
present to suppress H2 formation. The threshold n'^{}i2) is 
therefore higher than that of q = 1 for a given J21. This 
is the reason why the convergence of Xh2 is rather slow in 
Fig. In either case. Fig. ^ assures that the H2 phase 
is realized at the center of a contracting cloud before it is 
rotationally supported, as long as J21 is not far above the 
level J21 ~ 1 observed at z ~ 3. 

The above critical densities n'^{¥Ll) and n'^(H2) at 
J21 ^ 1 with a = 1 are in rough agreement with analytical 
estimations by Tajiri & Umemura (2000) for an isothermal 
sphere; n^{HT), n^{}i2) oc J2i^ ■ This confirms, in the case 
of a = 1, that nH(HI) is essentially the density at which the 
number of recombinations in a cloud exceeds that of ioniz- 
ing photons. Similarly, (H2) at J21 ^ 1 corresponds to the 
density at which the H2 cooling rate overcomes the UV heat- 
ing rate attenuated at the cloud center. For a much softer 
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Figure 7. Same as Fig.^, except for adopting a black-body spectrum with T^g 
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lO^K. In panels (a)-(c), all the lines overlap one 



spectrum (Tcff = 10*K), the deviations from such simplified 
estimations become larger, because of the greater signifi- 
cance of emitted photons and the LW photons. 



3.3 Criteria for collapse and H2 cooling 

We have performed the simulations for a variety of param- 
eters, such as the intensity and spectrum of the UVB, and 
the circular velocity and collapse redshift of a cloud. The 
results are summarized in the following. 

Fig. |l^ shows a J21 — K diagram on collapse and cool- 
ing of a cloud with Zco ~ 10 {zc = 3). The final state of the 
central gas is classified into three categories: (i) it collapses 
and cools to T < 5000K, i.e. the H2 phase is realized; (ii) 
it collapses but does not cool to T < 5000K; and (iii) it is 
prohibited to collapse within the present age of the universe. 
Also plotted for reference is the relation Tvir = Te™^", where 
T™''" is the maximum equilibrium temperature defined semi- 
analytically in Section 3.3 of Paper II. This relation roughly 
corresponds to the Jeans condition at the maximum expan- 



sion of the central gas and provides a reasonable estimation 
for the collapse threshold in the optically thin limit. 

The predictions of the full calculations are distinct from 
the optically thin results in two ways. One is that the col- 
lapse of low mass objects with Tvir < 10*K {Vc < 17 km 
s~^) is possible at small UV intensities. The other is that 
the central gas in collapsed clouds can cool efficiently to tem- 
peratures well below lO^K, as long as J21 is not far above the 
observed level J21 ~ 1. The suppression of H2 cooling occurs 
at J21 > 10" for a = 1 and J21 ^ 10^ for T^s = 10*K, mainly 
due to photodetachment of H~ as well as photodissociation 
of H2. To quantify this effect, we indicate in Fig. |l^ the av- 
erage intensity of the external UVB defined by equation 
at 0.74-13.6 eV. The suppression of H2 cooling is shown to 
take place at (J(0.74-13.6 eV)) 5 x 10* in massive clouds 
for both spectra. 

Given circular velocity and collapse redshift of a cloud, 
one can compute the critical radiation intensity for collapse 
and H2 cooling. Fig. ^ exhibits, as a function of the cen- 
tral collapse redshift Zco, the intensity above which H2 cool- 
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Figure 8. Same as Fig.^, except for showing a high mass collapse with Vc = 50 km s ^ in the UVB with J21 = lO'^ (long dashed), 1 
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ing is prevented completely in a cloud, i.e., the final state 
of the central gas belongs to the category (ii) or (iii) de- 
fined above. The critical intensity is in general an increas- 
ing function of redshift, i.e., gas density. It also increases 
with circular velocity but converges at 14 ^ 50 km s~^ and 
Vc > 30 km s"^ for a = 1 and T^ff = lO'^K, respectively. The 
convergence takes place once photoheating no longer affects 
the cloud evolution and only photodestruction of H2 and 
H~ becomes important. The converged intensities therefore 
have nearly identical average values below the Lyman limit 
(J(0.74-13.6 eV)) ~ 30{1 + Zcof for both spectra. Once con- 
verged, the critical intensities tend to drop slightly for larger 
Vc (panel b) because the electron fraction of the gas is lower 
in a cloud with higher hydrogen column densities. 

To explore the possibilities of different epochs of reion- 
ization, we have added to Fig. |l^ the results of runs with 
ztjv = 10 and 50. If Zco is sufficiently close to zuv, the 
central gas is already contracting when it is exposed to the 
UVB, and hence the stronger intensity is required to pre- 
vent collapse and cooling. The difference is greater for the 



harder spectrum and the smaller cloud mass, for which sup- 
pression of cooling is mainly via photoevaporation and de- 
pends highly on the initial dynamical state of the gas. For 
the softer spectrum and the larger cloud mass, on the other 
hand, the suppression is regulated to a greater extent by 
photodestruction of H2 and H~, and is less sensitive to the 
initial dynamical state. After all. Fig. |l^ indicates that there 
is a greater chance for the gas to cool even in a cloud well 
below the scale Vc = 30 km s~^ at high redshifts, provided 
that the intensity is comparable to or less than J21 ~ 1. 

Alternatively, given a specific model on the evolution 
of the UVB, one can deduce the critical circular veloc- 
ity or mass for collapse and cooling. Fig. ^ shows such 
quantities under which H2 cooUng is completely suppressed 
in an object with the cloud collapse epoch Zc- The UVB 
greatly prohibits cooling in clouds with Vc 30 km s~^ 
or A/cloud ~ 10^ M0 at low redshifts. The threshold values, 
however, depend highly on the time history and spectra of 
the UVB. For an evolutionary model given by equation (^, 
the critical mass drops by several orders of magnitude at 
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Figure 9. Evolution of (a) temperature, (b) Xh2 (thick lines) and Xhi (thin lines), versus gas density at the center of a collapsing cloud 
with Vc = 50 km s~^ and Zco = 10 {zc = 3) in the case of a = 1. Lines indicate J21 = 10^ (solid), 10 (dotted), 1 (short dashed), 10~^ 
(long dashed), and 10"'^ (dot dashed). For clarity, only the evolution at z < 15 is plotted. 



Zc ^1 3 (2cO ^ 10), enabling cooling in objects as small as 
M.ioud ~ 10*^ Mq. 

Fig. ^(b) also exhibits the mass and collapse red- 
shift of the 1,2,3(7 density perturbations in a cold 
dark matter (CDM) model with (fio, Ao, /iioo, ^b, fg) = 
(0.3, 0.7, 0.7, 0.05, 1.0). This is a specific example of the most 
successful cosmological model which is consistent with a 
number of current observations (e.g., Wang et al. 2000). Fig. 
|l4| shows that the UVB mainly suppresses collapse and cool- 
ing of ^ Icr objects with Mdoud ~ 10^ M© in the ACDM 
model. It is rather difficult to suppress cooling in higher-cr 
objects, unless the radiation intensity becomes much higher 
than the level observed at z 3. 

If different values of cosmological parameters are 
adopted, the collapse criterion in terms of circular veloc- 
ity remains practically unchanged as long as 57b /f2o <S 1. 
Thus, the critical mass shown in Fig. ^^b) scales with cos- 
mological parameters merely according to equation (^). 



4 DISCUSSION 

Our calculations predict that the formation of objects with 
Vc ^ 30 km s"'^ is greatly suppressed by the UVB at low red- 
shifts, whereas much smaller objects (Tvir ^ lO^K, T4 ^ 17 
km s-i) can collapse and cool at high redshifts {zc ^ 3, 
ZcO ^ 10). This is mainly because of three reasons: (i) in- 
tensity of the UVB tends to decline at high redshifts prior 
to complete reionization of the Universe; (ii) self-shielding 
of the gas is stronger at higher redshifts owing to higher gas 
density; and (iii) the gas density in Jeans unstable clouds 
rises rapidly along with run-away collapse, further promot- 
ing self-shielding and H2 formation. 

What are the fates of low mass objects that are able to 
cool and collapse in the presence of the UVB? If the cooled 
gas can promptly turn into stars, they are likely to constitute 
'building blocks' of present-day luminous structures. While 
such objects may be detected directly by future instruments 
such as NGST, they may also have links to observed sub- 
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galactic objects (e.g., Pascarelle et al. 1996) or dwarf galax- 
ies. In contrast, those prevented to collapse by the UVB will 
become gravitationally unbound and spread into intergalac- 
tic medium. They might contribute to Lya absorption lines 
in the QSO spectra (e.g., Umemura & Ikeuchi 1984). 

These speculations on the bifurcation of low mass ob- 
jects are illustrated schematically in Fig. uK The critical 
baryonic mass in the case of evolving J21 (eq!g) with a = 1 
is used to distinguish the uncoUapsed 'expanding clouds' 
from the others. Furthermore, we expect that the collapsed 
objects may be identified, mainly based on their distinct star 
formation histories, as at least two separate populations: the 
lower-mass ( ^ 10^ Mq), higher-cr ( Si Icr) systems and the 
higher-mass, lower-cr ones. The former population can be a 
plausible candidate of dwarf spheroidal/elliptical galaxies for 
the following reasons. First, star formation in these objects is 
expected to be active at redshifts higher than ~ 3, consistent 
with the fact that the majority of local group dwarfs show a 
sharp decline in the star formation rate at some 10 Gyr ago 
(Mateo 1998; Gnedin 2000a). Second, their gravitational po- 



tential is compatible with small velocity dispersions ( ^ 10 
km s"^) of local group dwarfs. Third, they are apt to be 
dark matter dominated and metal-poor since the gas infall 
and star formation are suppressed particularly at the outer 
envelope. Fourth, our simulations indicate that the central 
gas in these objects can cool efficiently by H2 in the course 
of spherical contraction before reaching the rotation barrier. 
Galaxy formation is thus likely to proceed rapidly in a dis- 
sipationless manner and lead to the formation of spheroids 
rather than disks. 

On the other hand, the latter population will exhibit 
signs of more recent star formation and contribute signifi- 
cantly to the global star formation rate of the Universe. The 
clustering of these low a objects is expected to be rather 
weak. They might thus have been identified as a part of 
faint blue objects (e.g., Tyson 1988; Cowie et al. 1988; Ellis 
1997) or star-forming dwarfs (e.g., Thuan & Martin 1981; 
Lee et al. 2000). This is in line with the idea proposed by 
Babul & Rees (1992) and extended by Babul & Ferguson 
(1996) and Kepner et al. (1997). One distinction is that the 
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impacts of the decline of the UV intensity at low redshifts 
turn out to be modest. This is because the collapse is de- 
layed by the kinetic energy of the gas attained during the 
peak of the UV intensity at 2: ~ 3 — 5 (although the thresh- 
old circular velocity does decrease with time, the threshold 
mass continues to rise due to a factor M/Vc ~ [1 + Zc]^'^^^). 

The present spherical collapse simulations are mainly 
applicable to low mass objects near the Jeans scale. Much 
larger objects, in contrast, are likely to first undergo sheet- 
like pancake collapse and evolve in a different manner from 
the spherical case. They would not end up with the run- 
away collapse as gravity of the sheet is readily overwhelmed 
by thermal pressure. In fact, Susa & Umemura (2000a,b) 
argued that the pancakes collapsing in the UVB can bifur- 
cate, depending upon the degree of self-shielding, into ob- 
served early and late type galaxies. Both their results and 
ours suggest that the UVB could play an essential role in 
regulating the star formation history of the Universe. 

Recent three-dimensional simulations show that a num- 
ber of subclumps tend to arise from the fragmentation pro- 
cess in a parent cloud (e.g., Bromm, Coppi & Larson 1999; 
Abel, Bryan & Norman 2000). The criteria obtained in the 
current calculations are applicable to each of such locally 
collapsing subclumps. As a net effect to the parent cloud, 
the gas dumpiness is likely to enhance self-shielding against 
the external UVB, since the number of photons necessary 
to ionize the gas increases roughly by a factor {n^) / {nn}'^ 
(Gnedin & Ostriker 1997; Ciardi et al. 2000b). 



Major uncertainties in our predictions are regarding the 
feedback from the sources formed in an collapsing object. 
While our simulations can probe directly the cloud evolu- 
tion until its central collapse, the evolution thereafter will 
be influenced by several additional processes. For instance, 
Omukai & Nishi (1999) argued that even a single OB star 
can photodissociate surrounding hydrogen molecules within 
a static cloud with TVir ~ lO^K. On the other hand, Fer- 
rara (1998) pointed out that blow-away processes during 
supernova explosions can enhance the fraction of hydrogen 
molecules in surrounding media. Such explosions may also 
disrupt low mass objects and induce metal injections into 
intergalactic space (e.g., Dekel & Silk 1986; Nishi & Susa 
1999; Ferrara & Tolstoy 2000). It still remains uncertain to 
what extent these processes affect successive gas infall and 
cooling in a dynamically contracting object whose central 
overdensity exceeds ~ 10^ . As this issue is beyond the scope 
of the present paper, we will investigate it separately in our 
future publications. 



5 CONCLUSIONS 

We have studied the formation of Population III objects ex- 
posed to the external UV radiation. Radiative transfer of 
photons, H2 formation and destruction processes, and hy- 
drodynamics are incorporated explicitly in our analyses and 
all shown to play essential roles in the development of Pop- 
ulation III objects. Although the UVB does suppress the 
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Figure 12. J21 — Vc diagram on collapse and cooling of a cloud with ZcO = 10 (zc = 3). Different panels compare the optically thin 
(left) and full (right) calculations, and the UV spectrum with q = 1 (top) and T^g = 10*K (bottom). Symbols indicate whether the 
cloud center collapses and the H2 phase is realized (open circles), it collapses but the H2 phase is not realized (open triangles), or it 
is prohibited to collapse (filled circles). Dashed lines show an analytical estimation for the collapse threshold in the optically thin limit 
defined in Paper II. The average intensity below the Lyman limit (J(0.74-13.6 eV)) is displayed on the right vertical axis. 



formation of low mass objects, the negative feedback turns 
out to be weaker than previously suggested in the following 
two respects. 

First, the cut-off scale of the collapse drops significantly 
below Tvir = 10*K {Vc — 17 km s"^) at weak UV intensities 
{J21 10^^), due to both self-shielding of the gas and H2 
cooling. Their importance becomes greater at higher red- 
shifts. At stronger intensities, objects as large as K ~ 40 
km s"i can be photoevaporated and prohibited to collapse, 
in agreement with previous investigations based on the op- 
tically thin approximation (e.g., Thoul & Weinberg 1996). 
High mass objects well above this scale can collapse almost 
irrespectively of the UVB by atomic cooling. 

Secondly, a spherical cloud that has become Jeans un- 
stable in the UVB tends to undergo run-away collapse and 
does not settle into hydrostatic equilibrium. As a result, 
even if the gas is ionized and/or H2 molecules are de- 



stroyed once by radiation at low densities (e.g., Haiman et al. 
2000), the radiation is attenuated and H2 can re-form dur- 
ing the dynamical collapse. The critical densities at which 
self-shielding and efficient H2 cooling take place depend sen- 
sitively on the spectrum and intensity of the external UVB. 
For a QSO-like {a = 1) spectrum, H2 formation is even en- 
hanced by a little radiation with J21 ~ 10""^ — 10"^ (Haiman 
et al. 1996b). For a stellar-like {Tc« = 10*K) spectrum, H2 
formation is hindered to a greater extent by the photons 
with energies below 13.6 eV (Haiman et al. 1997). In any 
case, the collapsing gas is shown to cool efficiently to tem- 
peratures well below 10'* K before rotationally supported, as 
long as the radiation intensity is comparable to or less than 
J21 ~ 1 observed at 2; ~ 3. The H2 fraction in the cooled 
gas reaches the level ~ 10""^, comparable to that achieved 
commonly in the metal-deficient post-shock layer (Shapiro 
& Kang 1987; Ferrara 1998; Susa et al. 1998). 
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Figure 13. Critical UV intensity above which H2 cooling is prevented as a function of central collapse redshift z^a'- (a) q = 1, and (b) 
Tcff = 10*K. Circular velocity of a cloud is taken to be Vc = 50 (thick lines), 30 (mildly thick lines) and 10 km s~^ (thin lines). The 
onset of the UVB is varied as zyv = 10 (circles, dotted lines), 20 (crosses, solid lines) and 50 (triangles, dashed lines). 



Our results imply that star formation can take place 
in low mass objects collapsing in the UVB. The thresh- 
old baryon mass for the complete suppression of collapse 
and star formation is Mdoud ~ 10^ M0 at Zc ~ 3 but drops 
significantly at higher redshifts. In a conventional ACDM 
universe, it coincides roughly with that of the la den- 
sity fluctuations at Zc ^ 3. Objects near and above this 
threshold can thus constitute building blocks of present- 
day luminous structures, and may also have links to dwarf 
spheroidal/elliptical galaxies and faint blue objects. These 
results indicate that the UVB can play a key role in regu- 
lating the star formation history of the Universe. 
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